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During  this  grant  period  ve  tackled  three  problems  all  related  to  syatema  with  rugged 
snergy  landaeapet  (1)  Developing  replica  molecular  dynamlsB  method  to  detect  bottle'- 
steks  in  energy  aurface;  (2)  Characterising  the  aelf^erganlsatlbn  of  biomoleeules  in 
Kerms  of  sequences;  (3)  Taming  the  energy  landacape  for  design  of  bionaterlals. 
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Proposal  Titls;  Explorin«  Rugged  Energy  Landscape  in  Large  Systems 
Pinal  Report 

Spoa»r.  Ail  Fteee  OffiM  of  Sdtttific  R«««k 

Grant  Nuznber:F4963O041OlO6 

PriMlp.1  In«.tig.tot:  D.  Thimmriol.  U»iWy  of  M«yl^d,  CoUog. 

MD  20742 

Th.  |««d  *° 

™tbod.  .0  «b«o=.«i»  a-  -re.  and  «ddk 

p<dats)  ia  tatSf  "*  “ 

biotojy.  ««1  «i«=.  in  wMd.  on.  i.  int^tod  in  the  dy«». 

ics,  »»t«W  «.d  of  .,.(«» in  hieh  dtoondooe.  E  iampta.  whae  thi. 

the  no®  indnd.  kineti-  of  fbiang  of  ttonoleeuk..  djnnunc.  <d  the  li,oid  to 

^  ttandtlan,  behsviot  of  poljmen  et  inUifccM.  fat  ail  hoee  ptoblems 
.b««tic.l  de«iption  of  the  v«i«.  p™?«ti-  «“  datemdnetion  of  the 

t..„sitioert.*..inhi^din»ndonalfe«  «.«».«*«•  Dotitt  the  px««nt  fond- 
ing  p«iod  w  «complished  three  import«t  p»l.  rdeted  t<.  out  undetmeadin* 
uee  of  es-rteme  whoe.  und«lyin*  energy  londteepe  »  cotnjd-t.  Throu^out 
<„„«..rcl.«u«d««.Vl..fcomd»kin.ti=.ofproteinf<^^  Themethode 
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developed  are  quite  geaeral  and  should  be  applicable  m  tbs  d«^  of  novel 


(a)  Development 
Ibarner  hdlghts  in 


polym«ie.ubrtfflc«. 

rf„pUc»  »ol«»kr  dyMiide.  (RMD)  to  ptob,  the  ainia.  »nd 
systems  Kith  »  distribution  o£  minimi  end  msjdme-  All  of  the  sxemples  that  we 

h.v,c«»id«.dbeiongt.lil.»t,ory.  (b)Theu..o{.n*rgy1»d.cpeid...to 
ds«ifyfcldingkin.ti«cfpiot.in..  (cjDsrign  of  peptide  «,««««.  that  Wd  into 

e  J-sifod  staueluie.  A  brief  description  of  each  fellows. 

PTr«-  Molecular  Dynamice:  Consider  a  rough  energr  todscspe.  Here 

•  eze  eeoarated  from  eaeb  other  hj  barriers  that  are 

there  are  many  mmima  that  are  separatea  mva** 

loea«d  only  on  long  time  ecti...  I.  typical  ti-ulath®*  «f 

that  the  computational  time  seal,  ad^prat.  .«n,  Jing  of  the  emdor- 

snatimw  take  pUce.  ThU  U  the  etgodic  theorem.  If  the  rimuUtlon  tmnperatui. 

is  not  high  «»mgh  fer  the  barriem  to  he  overcome  tl»n  ti»  .ystem  ^ 

trapped  in  one  of  the  rtunim.  and  relUbl.  results  cannot  be  obtrined.  The  first 

question  thet  ...d.  to  be  enmred  is  hms  th.  mdst-se.  of  tie  multir^^ 

tore  be  detected?  W.  de-«lop«i  a  novti  rimulation  mathoc,  trfeirwi  to  a.  the 
replica  molecular  dynemim  (HMD),  in  which  «o  replicu  of  he  system  with  rim 
distinct  initial  condition.  generrted.  For  each  initial  cotditice.  «.  calculate 
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that 


ivsrsa 


the  time  average  energy  of  izuiividual  particles  assmniag  that  thej 
ftmetion  of  th©  system  is  given.  If  the  observation  time  sale  is 
value  is  the  average  energy  of  the  given  particle  wbidi  is  a  vel, 
dynamic  quantity.  We  introduced  the  energy  metric  which  ea 
dependence  of  the  distance  between  the  time  averaged  values  of 
particles  starting  with  the  two  distinct  replicas.  We  showed 
which  the  confonnational  space  is  sampled  efficiently  the  in 
metric  goes  as  Dt,  where  D  is  the  rate  of  eqiloration  ci  co 
there  is  a  bottleneck  between  the  two  basins  of  attraction  then 
eturatee.  The  temperture  dependence  of  the  energy  metric  giv! 
xier  hei^  separating  the  basins  of  attraction.  We  have  appUec 
vndexstand  the  dynamics  of  proteins.  In  particular  we  above  i 
there  are  and  barriers  rangibg  from  the  small  energy  s 

scales.  Thus  in  order  to  understand  macromolecular  sel^cag 
account  for  motions  from  few  angstroms  length  scales  to  th« 


infbixaa' 


system. 


des 


Folding  Kinetics  of  Proteins:  A  concrete  test  of  these  i 
the  folding  kinetics  of  proteins.  Proteins  are  known  to  be  gov^i 


potential  energy 
the  long  time 
defined  thenso* 


.culated  the  time 
the  energy  of  the 
for  syst^  in 
of  the  energy 
.tional  space.  If 


fthe  energy  metric 
I  the  average  bar- 
these  methods  to 
rigouroua^  that 
4aleB  to  the  largest 
laatisatiai  one  has  to 
size  of  the  entire 


is  to  understand 
ed  by  a  complex 
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energy  laadacape.  In  ot<Ur  to  under#  W  the  kinetics  in  terms 
tou^  energy  surface  and  connect  to  experiments  it  is  nectts 
properties  of  the  landscape  to  intrinsic  observables  of  piuteinj 
that  the  folding  characteristic*  can  be  sucdatly  accounted  for  ty 
teristic  temperatures  One  of  them  is  T*  which  is  the  temperat 
ia  a  confonaational  collapse  feom  an  extended  confonnatlon 
The  other  is  at  which  the  protan  adopts  the  unique  nativej 
that  minimisation  of  the  difference  in  these  two  temperatures 
landscape  thus  guiding  the  pcdypeptide  chain  into  the  uni<iue 
vBTsely  if  the  difference  is  large  then  the  landscape  is  rough  wbi 
folding.  The  major  «mclusion  is  that  ail  the  complexities  of 
ergy  surface  can  be  folded  into  this  one  sln^  parameter  and 
folding  kinerics  emerge  by  merely  tuning  this  difference  in  t 
Design  of  folding  Sequences:  If  it  is  correct  that 
two  characteristic  temperatures  leads  to  fast  folding  in  a 
surface  then  it  must  be  possible  to  design  optimiaed  seqaetj( 
into  a  given  structure  i.e.,  one  should  be  able  to  apprccdmai 
folding  problem.  We  set  out  to  do  this  m  the  contact  of 


pf  the  underlying 
to  connect  the 
We  discovered 


the  two  charae> 

:e  at  which  there 
a  compact  state, 
state.  We  showed 
l^ds  to  a  smoother 
native  state.  CSon* 
Pleads  to  slugpsh 

the  tou^  free  ea* 
'^larious  scenarios  for 
iBihperatures. 

X  ainimization  of  the 
ir  Mother  free  energy 
ces  that  would  fold 
.ttiy  solve  the  inverse 
design,  of  four  helix 
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buadlw.  It  hM  been  ihownexperimeatallythit  Urate  ' 

rf  andnoraid.  that  trould  fidd  into  fcui  babx  bmidl«.  We  rtatj 

,nadd  contUnta*  jn.t  thm  rarfnoiAl  teddtt- (Mo«^ 
w,  petformwi  matte  cjdo  in  reqnenee  epece  end  nramtoted 

yielded low-valoetrf the (fifiaence between T. end Tp.  The 

wore  a  lew  eequencee  that  ean  yield  opUalted  etruetoiee.  We 
Idartiea  amnintora  to  riiow  that  theie  eequencee  indeed  feW 

taming  the  i<i«ed  landee^te  we  beee  come  up  with  a  dedgn 
»qoeaeee  that  fold  into  a  ^  itturture.  We  bdieve  that  tbii 

and  applienUo  {or  dedgun*  pdymerie  oateriele  on  meeoeeopje 


.did&te  sequesces 
Ud  ^th  A  simple 
).  THen 

sequeaicee  ihftt 


exaeriinexsts) 


,t  was  thsit  there 
then  pez£ornied  a 
rapidly.  Thus  by 
iterion  for  finding 
method  is  general 
length  scales.^ 
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